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^ ; Abstract 



We investigate the relation between non-unitarity of the leptonic mixing matrix and leptogenesis. 



We discuss how all parameters of the canonical type-I seesaw mechanism can, in principle, be 
q^, reconstructed from the neutrino mass matrix and the deviation of the effective low-energy leptonic 

r-j | mixing matrix from unitary. When the mass M of the lightest right-handed neutrino is much 



lighter than the masses of the others, we show that its decay asymmetries within flavour-dependent 
leptogenesis can be expressed in terms of two contributions, one depending on the unique dimension 
five (d = 5) operator generating neutrino masses and one depending on the dimension six (d = 
6) operator associated with non-unitarity. In low-energy seesaw scenarios where small lepton 
number violation explains the smallness of neutrino masses, the lepton number conserving d = 6 
operator contribution generically dominates over the d = 5 operator contribution which results in 
a strong enhancement of the flavour-dependent decay asymmetries without any resonance effects. 
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To calculate the produced final baryon asymmetry, the flavour equilibration effects directly related 
to non-unitarity have to be taken into account. In a simple realization of this non-unitarity driven 
leptogenesis, the lower bound on M 1 is found to be about 10 8 GeV at the onset of the strong washout 
regime, more than one order of magnitude below the bound in "standard" thermal leptogenesis. 
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I. INTRODUCTION 



Non-unitarity of the leptonic mixing matrix at low energies is a generic manifestation of 
new physics in the lepton sector often related to the mechanism responsible for the generation 
of neutrino masses. Non-unitarity appears whenever additional heavy particles mix with the 
light neutrinos or their charged lepton partners. After integrating the heavy states out of 
the theory, the 3x3 submatrix of the light neutrinos remains as an effective mixing matrix. 
This low-energy leptonic mixing matrix is, in general, not unitary. 

One example where non-unitarity is predicted is the generic type-I seesaw mecha- 
nism , where the Standard Model (SM) is extended by (typically three) right- 
handed neutrinos. If the type-I seesaw mechanism operates at energies as high as the Grand 
Unification scale (GUT scale), then non-unitarity effects are tiny. However, if the seesaw 
mechanism is realized at low energies close to the electroweak scale, then non-unitarity is 
enhanced and can be observable. It may then provide important hints to the origin of neu- 
trino masses. It is important to note that, in low-energy seesaw scenarios, the smallness of 
the neutrino masses is not explained by the largeness of the seesaw scale. However it can be 
explained in a technically natural way by a lepton number symmetry which is broken only 
by a small amount Q]. 

While the neutrino masses in the type-I seesaw mechanism are effectively described by the 
unique lepton number violating Weinberg operator of d = 5, the non-unitarity of the leptonic 
mixing matrix is generated by the lepton number conserving d = 6 operator contributing to 
the kinetic terms of the neutrinos . In low-energy seesaw scenarios with approximately 
conserved lepton number, the d = 6 operator can cause significant effects, since it is not 
suppressed by the smallness of the neutrino masses. 

One attractive feature of the seesaw mechanism is that it can explain the observed baryon 

, 1 

asymmetry of the Universe via the mechanism of leptogenesis [8j (for a recent review, see 
Ref. [9|). If the seesaw mechanism, and thus also the mechanism of thermal leptogenesis, 
operates at high energies, the decay asymmetries for leptogenesis are typically dominated by 
the d = 5 operator. Unfortunately the high-energy parameters which control leptogenesis 
cannot be fully reconstructed from the measurements at low energy, since combinations of 
the Yukawa couplings different from those in the d — 5 operator also appear. 

On the other hand, if the seesaw mechanism operates at lower energies, predicting an 
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observable non-unitarity of the leptonic mixing matrix, one may in principle obtain enough 
information from low-energy measurements to reconstruct the full Lagrangian and, therefore, 
the parameters that control leptogenesis [ill. 



The purpose of this paper is to clarify the relation between the high-energy parameters 
that control successful leptogenesis and their low-energy manifestations, i.e., neutrino masses 
and mixings and deviations from unitarity of the leptonic mixing matrix, making special 
emphasis on the latter. In Section [III we describe in detail the method and conditions under 
which the full high-energy Lagrangian can be reconstructed from the low-energy effects. In 
Section II I II we introduce a "minimalistic" low-scale seesaw model with three right-handed 
neutrinos in which the smallness of neutrino masses is explained by an approximate lepton 
number symmetry. However, deviations from unitary mixing induced by the d = 6 operator 
are not protected by the symmetry and can be sizeable, leading to effects that could, in 
principle, be tested in precision electroweak measurements. In Section [IV] we discuss how, 
in this model, leptogenesis could be driven by the d — 6 operator that induces the deviation 
from unitary mixing via flavoured leptogenesis. In Section [V] we point out that the same 
d = 6 operator that drives flavoured leptogenesis can also lead to a flavour equilibration, 
which could wash out the generation of lepton number. We also discuss the conditions under 
which successful leptogenesis can occur. Finally, in Section IVIl we summarize and discuss 
our results. 



II. NON-UNITARITY RELATION TO HIGH-ENERGY OBSERVABLES IN THE 
TYPE-I SEESAW MODEL 

In this section we describe how the full Lagrangian of the type-I seesaw can, in principle, 
be reconstructed from the low-energy observations of neutrino masses and mixings, including 
deviations from unitary mixing. The conditions under which this reconstruction is possible 
were described in Ref. J3], while a method to realize the reconstruction was outlined in 



Ref . [lOj . Here we present a new algorithm to perform the reconstruction through which the 
high-energy parameters can be derived more easily. 

Let us consider the Lagrangian of the standard type-I seesaw model which consists of the 
one for the SM plus an extra piece containing the allowed couplings between the SM fields 
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and additional gauge singlet fermions (i.e., right-handed neutrinos) N^: 

J? = J?sm~ ^NiM^N c i - (Y N ) ia Ng>Hl + H.c. . (1) 

Here, <fi denotes the SM Higgs field, which breaks the electroweak (EW) symmetry after 
acquiring its vacuum expectation value (vev) t>Ew> an d we have used the definition = ir 2 0*. 

The low-energy effects of the three-family low-scale seesaw model, from the point of view 
of neutrino oscillation experiments, is given by two effective operators, one of mass dimension 
five and one of mass dimension six. The d = 5 operator is the ubiquitous lepton number 
violating Weinberg operator 



6C d = 5 = - c d a f (lcJ*) [9 LA + H.c. , (2) 



1 ~. 
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which is the lowest- dimensional effective operator for generating neutrino masses using the 
field content of the SM. The coefficient matrix c^ 5 is 

= -(Y^) ai (M N )^(Y N )^ (3) 

and relates to the low-energy neutrino mass matrix as 

m v = vl w c d=5 . (4) 

The effective d — 6 operator 

5C d=& = c% & (Zjy0^L^ (5) 

conserves lepton number and, after EW symmetry breaking, contributes to the kinetic terms 
of the neutrinos. After their canonical normalization, they generate a non-unitary leptonic 
mixing matrix iV as well as non-universal couplings proportional to N*N of the neutrinos 
to the Z boson (see, e.g., Refs. js, 7, 12]). The coefficient matrix c^ 6 is given by (see, e.g., 



Ref. fl) 

C S 6 = J2( Y NU( M N)u 2 (Y N ) ip , (6) 

i 

in the basis where Mjv is diagonal. 

If we parametrize the non-unitary leptonic mixing matrix iV as [3] 

N=(1 + V )U, (7) 



where rj is Hermitian and U is unitary, then r] a p is related to the coefficient matrix c^ 6 by 

Vafi = ~v 2 EW cif/2 . (8) 

In Refs. I12I, 14, 15, Q, the following constraints on these parameters at the 90 % C.L. 
were derived: r] ee < 2.0- 1(T 3 , r/ eM < 5.9- 10" 5 , r] eT < 1.6 -lO" 3 , ij^ < 8.2- 1(T 4 , rj^ < 1.0- 10" 3 
and 7] TT < 2.6 • 1(T 3 . 

In the type-I seesaw the full 6x6 mixing matrix U tot is the unitary matrix that diagonalizes 
the extended neutrino mass matrix: 
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U tot = | 


' m 













M N J 




V 


M ) 



UL U t/tot = , (9) 

\m D M N J \ M J 

where = vew^n and are the neutrino's Dirac and Majorana mass matrices and m 
and M are diagonal matrices. It is easier to perform the diagonalization in two steps: first 
a block- diagonalization and then two unitary rotations to diagonalize the mass matrices of 
the light and heavy neutrinos, i.e., 

(a b \ ( u \ 

u«* = , (10) 

\cd) \o vj 

where U and V are unitary matrices. Without loss of generality, we can choose a basis for the 
heavy singlets such that V — I and Mjy = M. When performing the block diagonalization, 
the mixing between the light and heavy neutrinos is suppressed when > mp so that 

B ~ 6 = m ] D M N l . (11) 

We can exploit the suppression of Eq. ffTTT) to write the unitary block diagonalization as 
the exponential expansion of an anti-Hermitian matrix: 

( e\ ( i-ieet e \ 

= exp = 2 +0(6 3 ). (12) 

y-et J y -et 1 - |ete y 

The block-diagonalization yields the complex symmetric neutrino mass matrix 




-m^M-^o, (13) 



which can be diagonalized by a unitary transformation U such that m = diag(m 1? m 2 , m 3 ) = 
U T m v U . Notice that the mixing matrix of the three light neutrinos is given by 

N = AU = (1 + ij)U = (1 - &&y2)U, (14) 
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As described above, 77 = —00^/2 exactly contains the coefficients c d=6 of the d = 6 operator. 
In particular, note that this implies that c d=6 is Hermitian and positive semidefmite. 

Assuming that the low-energy observables U, m and 77 have been measured, it is then 
natural to ask the question of whether one can reconstruct the high-energy parameters 
contained in and M. In order to split the low- and high-energy observables, we make 
use of the parametrization proposed in Ref. [18( for the Dirac mass matrix mo, which is 



introduced as follows: From Eq. ffl3|) we have that — \J m x U T m r D M x muU\' 



m 



-1 



Thus, defining R = iy M~ l m£,U V m _1 , we have the condition R T R = 1. We will consider 
here the case in which the number of heavy right-handed singlets is equal to the number of 
light neutrinos, and thus, the matrix R must be a complex orthogonal matrix. Multiplying 
by inverses of matrices, the definition of R can be rewritten as 



m D = -i^MRy^RU ] (15) 



Notice that the matrix U that appears in the parametrization of Eq. (1151) is not the neutrino 
mixing matrix that describes the neutrino couplings in charged-current (CC) interactions 
N = (1+7])U, but only the unitary part of the CC mixing diagonalizing m u . In Eq. ffl5l) the 
observables of the d = 5 operator describing the light neutrino mass matrix are all contained 
within U and m, while R and M contain the missing information in order to reconstruct 
the high-energy parameters. In order to perform the reconstruction we will assume that 
both the d — 5 operator, i.e., U and m, and the d = 6 operator rj are known. Notice that 
the deviations from unitary mixing encoded in rj can be probed and potentially measured 

as well as neutrino oscillation experiments 
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ectroweak decays [12J, 
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Using the parametrization of Eq. (fl5l) in the expression for the d = 6 operator 2r] = 
—m} D M~ 2 niD, we define the matrix H as 

H = -Vr^U^r/UVr^ = {R*y 1 M~ 1 R. (16) 

With this definition, if is a positive semidefinite Hermitian matrix, since it can be decom- 
posed as FF\ and contains all of the available low-energy information. Equation ( fl6l) is 
known as the conjugate diagonalization (or simply the "condiagonalization" ) of the matrix 
H. Notice that replacing (R*)" 1 by R' 1 in Eq. (fl6|) would reduce it to a normal diagonaliza- 
tion. It can be shown that all Hermitian positive definite matrices can be condiagonalized 



and the solution is unique under the requirement that the M{ are real and positive. The 
simplest way of reconstructing R and M from the matrix H is to note that H*H = H T H is 
a complex symmetric matrix and that 

H*H = R~ 1 M~ 1 R*(R*)~ 1 M~ 1 R = R T M~ 2 R, (17) 

and thus, R is the complex orthogonal matrix which diagonalizes H*H with corresponding 
eigenvalues M~ 2 . It should be noted that the reconstruction of R and M has been previously 
studied in Ref. [h]]. However, the simple reconstruction algorithm presented here, only 
involving the diagonalization of a complex symmetric matrix, is new. Notice that, for the 
algorithm to work and for Eq. (Tl6]) to have solutions, it is vital that the matrix if is a 
positive definite matrix. This will be guaranteed as long as the d = 6 operator 

c d=e ig 

also 

positive definite. By construction positive semidefinite so the only case in which H 

is not positive definite is when c d=& has one or more zero eigenvalues. This could happen for 
two reasons: either the number of heavy right-handed neutrinos is smaller than the number 
of light neutrinos, or the rows of the Yukawa matrix are not linearly independent. In any 
of these cases the number of parameters in the high-energy theory is smaller than in the 
full case considered here and the reconstruction of the Lagrangian would then be easier, 
sometimes without the need of involving the d = 6 operator but only through U and m (see, 



e.g., Ref. |25|])- If, on the other hand, the number of heavy right-handed singlets is larger 
than the number of light neutrinos, then the information encoded in the d — 6 operator 
7] plus the neutrino masses and mixings m and U is not sufficient to reconstruct the full 



high-energy Lagrangian. We find these limitations to be in agreement with Ref. 
the reconstruction algorithm relies on rj being invertible. 



3, 



where 



III. THE THREE-FAMILY LOW SCALE SEESAW SCENARIO 

The "minimalistic" low-scale seesaw model we present here is a type-I seesaw model 
with three right-handed neutrinos to which we additionally impose a softly broken "lepton 
number"-like (global) U(l) symmetry where the charge of the SU(2)l doublets is opposite 
to that of the field but equal to that of the field N^. We assign zero lepton number to 
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./Vr. In the symmetry limit, and Y/v are forced to have the form: 



^0 0^ 



N 



Ve Vn Vt 





/ 



M 



N 



M' 
M 
M 



(18) 



At this level, the neutrinos are exactly massless, but non-unitarity of the leptonic mixing 
matrix is already induced. When small soft breaking terms fii and \J i are allowed, this rigid 
structure is perturbed to 



Y, 



N 



Ve P'n 
Ve Vfi Vt 



M 



N 



/4 Pi M 
\p' 5 M fi 5 J 



(19) 



and masses for the light neutrinos, suppressed by the small pi and p! i parameters, are gener- 
ated. Notice that the d = 6 operator is not protected by the lepton number symmetry and 
would have a large leading order contribution 

VaV/3 



d=6 



+ ..., 



(20) 



a/3 ' M 2 

where the dots denote sub- leading terms proportional to p a (or p^/M, p 5 /M) and p' a (or 
p' 4 /M, p' 5 /M). In particular, the six moduli \r} a p\ depend only on the three parameters 
\y e /M\, \y fl /M\, \y T /M\. This implies that at leading order the rank of the d = 6 operator is 
one and that the reconstruction algorithm presented in Section [III would not be applicable. 
On the other hand, the number of parameters contributing to the d = 5 and d = 6 operator 
at leading order in this model is very limited and they can actually be reconstructed with 



information mainly from m and U . Indeed, the d = 5 operator is given by 



25): 



d=5 I ^5x1 . 1 , Pb x 



(21) 



Following Ref. [25] both vectors y a and p a — {Hb/M)y a can be reconstructed from m and U 



up to an overall normalization. In particular, 



y a oc sj\ + pf/* 3 + 
P 



pu ( 



all 



\/T+r - y/r 



a/1 + r + y/r 



(22) 
(23) 



for normal hierarchy and 



oc VTT~pU* a2 + y/T^~pU* al , 



(24) 




FIG. 1: Effective operator decomposition of the diagrams leading to leptogenesis. 



y/T+r-1 
VI + r + 1 

for inverted hierarchy. Here, r = | Am^l/I Amjy. If terms beyond leading order are consid- 
ered for the d = 6 operator, then the reconstruction described in Sec. [Ill can be applied. In 
principle, this would also allow the extraction of the first Yukawa row Y al = /jf a if the d — 6 
operator is known with sufficient precision. 



IV. LEPTOGENESIS IN THE THREE-FAMILY LOW-SCALE SEESAW SCE- 
NARIO 

A. Flavour-dependent decay asymmetries 



The flavour-dependent CP asymmetries for the lightest right-handed neutrino (with mass 
M') in the lepton flavour a are given by 



87r(YYt) u 



E Im 



(26) 



where Xj = Mf/M* > 1 was assumed (with j = 2, 3, M x = M', and with M 2 , M 3 being the 
masses of and iV|J. It is then possible to rewrite S\ i0L as the sum of two contributions, 
that of the d — 5 operator and that of the d = 6 operator (see Fig. [1]): 



8vr(YYt 



Hi 



3M' 



y* „d=5 yt _|_ liyf'jS y d=b y 



12 - 



*/9 1 /31 



M' 



/2 



87r(yyt) u 



E{ Im [^ c S 6 4]} 



(27) 



where, in the last step, we have neglected the contribution from the d — 5 operator, since it 
is protected by the lepton number symmetry that explains the smallness of neutrino masses, 



which is not the case for the d = 6 operator. Assuming the presence of an order one phase 



|Y"i T |, we can further 



contribution to the imaginary part in Eq. (J27J), as well as |Yi e | 
simplify this expression and express it in terms of the non-unitarity parameters in Eq. 

M' 2 x ^ d=6 1 M' 2 



£1, 



E 



24tt ^ ^ P 

P^a 



12n v 



EW 



p^a 



(28) 



For M' = 1 TeV and r] a p of order 10 4 we obtain a large CP asymmetry £ 1)Q of order 
10 -3 without resorting to the usual enhancement for quasi-degenerate masses of the heavy 



neutrinos 
Ref. 
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27| . As the contribution to £i )Q originates from the d = 6 operator (see also 



28|). we will refer to this scenario as non-unitarity driven leptogenesis. 



It is crucial for our scenario that leptogenesis occurs when flavour effects are rele- 
vant 29|, |30j, i.e., for a right-handed neutrino mass M' < 10 12 GeV, so that the flavoured 
CP asymmetries need to be considered instead of the total one S\ = ^2 a £i, a , which is sup- 
pressed by the neutrino mass operator (d — 5). As we will see next, another requirement for 
a non-zero baryon asymmetry is that the washout in each flavour is different. This is there- 



fore an example of purely flavoured leptogenesis 
exclusively due to flavour effects. 



3d, 



3JJ, where the asymmetry is generated 



B. Final baryon asymmetry 

An important parameter for flavoured leptogenesis is given by the decay (or washout) 
parameter, Ki a , induced by the RH neutrino Ni, defined as the decay width over the Hubble 
expansion rate when T = Mf. 



r D (A^40 + 40t) = \ Yia \ 2 v 2 



K __ J- ZH iv » ~^ ^aV "I" ) _ \ I ia\ ^EW (OQ") 



where m* ~ 1.08 x 10~ 3 eV [32|. For future use, we also define Ki = J2 a ^ia- 

With the assumption M' -C M and M' < 10 9 GeV, the asymmetry is generated by Ni in 
the three-flavour regime The contributions from the heavier states are exponentially 

washed out in all flavours, as long as Ki a > 3, Vet G {e,yU,r} |33f ] . Moreover, as pointed 
out above, since the CP asymmetry is non-zero due to a pure flavour effect, it is crucial for 
the generation of the asymmetry that the washout is different in each flavour. Therefore, 
for simplicity, we will assume that the asymmetry is dominantly generated in one flavour, 
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i.e., r/Ba — Vb- The flavoured final asymmetry can be expressed as 



r] Ba = 0.88 x 10 2 e ha n(K la ), 

where it was assumed that sphalerons decouple after the electroweak phase transition 
For Ki a > 3, the efficiency factor k can be approximated as [36] 

0.5 



(30) 



34 



35] 



la ) 



1.16 ' 

a 



(31) 



In order to focus on the more relevant non-unitarity parameters and the mass scale M', we 
will for now fix a washout K\ a = 5, which is a typical value in the strong washout regime. 
This requires Yukawa couplings \Yi a \ ~ 10 -6 , as needed for the TeV-scale seesaw mechanism. 
We further assume the value K\ = ^ 7 -^i 7 = 45, such that K la <C K\^ a . In other words, 
in terms of Yukawa couplings, we have the relation 2 |Yi a | = (Yi^a]. Using Eqs. (1271) . (1301 



and fl3TT). we have 



0.88 x 10- 2 M' 2 ^ ,n- 5 f M ' 

2_^Vaf3 - 0.6 X 10 I — 



367T 5 



1.16 



M 



VaVp- 



(32) 



This prediction should be compared to the measured value r]g MB = (6.2 ±0.15) x 10~ 10 |3 
With Yukawa couplings y@ ~ Oil), it is easy to see that leptogenesis is possible with a mild 
hierarchy M'/M ~ 10~ 2 . Moreover, it is interesting to see that the scale of leptogenesis can 
be lowered, at least in principle, to the weak scale. This is possible in our scenario because the 
purely flavoured contribution to the CP asymmetry is not suppressed by the neutrino mass 
operator as in the usual case. The lower bound on the scale of leptogenesis 
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39, 



40 



411, 



given by Mi > 3 x 10 9 GeV at the onset of the strong washout, therefore does not apply. 



This was already noticed in Ref. 



331 ]. where, using the Casas-Ibarra parametrization 



18| 



in the limit ^32] 3> 1, it was shown that the scale of leptogenesis could be lowered; the 
inverse seesaw model can be shown to correspond to the extreme case ja^l ^> 100. The 
seesaw model under consideration could therefore potentially offer an al 
degenerate RH neutrinos 0, 27] to evade the gravitino bounds 
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43 



ternative to quasi- 
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V. FLAVOUR EQUILIBRATION 



From the previous discussion it seems that the scale of leptogenesis could be lowered to 
the weak scale without any problem. However, an important effect was neglected, namely 

11 



flavour equilibration [46J. If flavours equilibrate, the final baryon asymmetry is proportional 
to the total CP asymmetry e\ = J2a £ ^a, which is suppressed by the d = 5 neutrino mass 
operator. We would thus recover the standard scenario, and the usual lower bound would 
apply 

Let us now estimate how efficient flavour equilibration is in our case. The main processes 
are AL = scatterings with off-shell N2 and N3, e.g., £ a <f) — > tp§. Contrary to standard 
leptogenesis, the rates can be large, since the Yukawa couplings y a are not constrained by 
neutrino masses, and therefore flavour equilibration is potentially a problem. The question 
is to what extent it reduces the available parameter space. 

There are three different channels contributing to AL = scatterings: s-channel £ a <p — > 
), t-channel £ a <ft — > iptft, and t-channel £ a £% — > (fxp^. The reduced cross-sections for these 



47J], and in the limit M' -C M the total AL = cross-section 



Mi) ^wwm:v (33) 



processes can be found in Ref. 
is given by 

4 7T \M J 

where x = s/M' 2 . Note that the N 2 and N 3 contributions are essentially equal. The reaction 
rate is then obtained using 

m'z 2 r°° 

r*t° = 967r2c(3) J dx^izVx)^), (34) 

where z = M'/T and /Ci is the modified Bessel function of the second kind. This rate, as 
all the rates entering the Boltzmann equations, will be compared to the Hubble expansion 
rate, given by 

M' 2 

H(z) = IM^g;——, (35) 
z z Mp\ 

where g* = 106.75 and M P1 = 1.22 x 10 19 GeV. 

The Boltzmann equations for leptogenesis are given by 
dN Nl 



dz 



D(N Nl -N e £), (36) 



^ = E la D (N Nl - N%) - N Aa - ^ S^=°(N Aa -N Af} ), (37) 

where Nx denotes any particle number or asymmetry X in a portion of comoving volume 
containing one heavy neutrino in ultra-relativistic thermal equilibrium, so that N^(T ^> 
Mi) = 1. As a function of z, the equilibrium RH neutrino number density is given by 
N]f = 0.5 z 2 K,2{z). Furthermore, we have defined A a = AB/3 — AL a . The normalized 
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decay rate is given by D = Td/[H{z) z] = K\z (I/7) with the thermally averaged dilation 
factor (I/7) given by the ratio of the modified Bessel functions K\{z)lK,i{z). Finally, we 
defined the normalized inverse decay rate W^P = Y l ^/[H{z)z\ = 0.25 K la K\(z) z 3 , and 

-<AL=0 — -pAL- 
>a,p — 1 a/3 



S^r = °/[H(. z ) A- Note that the normalized scattering rate is fitted within 10 % by 



S*t° - 6.5 x 10- (^) z~\ (38) 

In the Boltzmann equations above, spectator processes Q, 13] and the conversion of a 
lepton flavour asymmetry into a A a asymmetry have been neglected, but we have checked 
that they do not change our results by more than 20 %. We have also checked that AL = 1 
scatterings contribute subdominantly in the strong washout regime under consideration. As 
for AL = 2 processes their rates are suppressed by the small neutrino masses and can be 
safely neglected. 

We show in Fig. [2] how the scattering rate normalized to the Hubble expansion rate 
varies with M', for a fixed hierarchy M/M' = 10 and washout K la = 3.5 , and for Yukawa 
couplings such that the baryon asymmetry of the Universe is produced [c.f. Eq. (|32|) ]. The 
shaded regions denote the asymmetry production time, roughly when zb — 2 < z < zb + 



2 with z B {K) ~ 2 + 4K 0A3 exp(-2.5/K) ~ 5 5_0J. Note that even though the flavour 
equilibrating scattering rate falls out of equilibrium before entering the shaded region for 
M' > 3 x 10 8 GeV, there is still a residual effect, which suppresses the asymmetry by a 
factor of 2-3. The reason is that the scattering rate is within a factor of 2 of the inverse 
decay rate during all the asymmetry production, and therefore some flavour equilibration is 
achieved. 

Since T AL=0 /H oc 1/M' for a fixed mass hierarchy M'/M, it is clear that the scale of 
leptogenesis cannot be arbitrarily low. We have solved numerically the system of Boltzmann 
equations (|36j) and (1371) without neglecting the contribution from any flavour, so that 

7] B = 0.88 X 10~ 2 ^ £ l,a K(Kla), (39) 

a 

and we have required that the final asymmetry should fall within the 3a range of the observed 
baryon asymmetry, i. e.,r]B > 5.75 x 10~ 10 . We then calculated the lower bound for successful 
leptogenesis for different hierarchies M/M', and the result is shown in Fig. [3J We find that 
the lowest bound on the scale of leptogenesis is given by 

MLn > 10 8 GeV, (40) 

13 



100 



M' = 


3xl0 7 


GeV 


M' = 


3xl0 8 


GeV 


M' = 


3xl0 9 


GeV 


.... W ID 

a 







&3 




Z = MIT 

FIG. 2: Scattering rate S^ =0 leading to flavour equilibration vs. z = M'/T for different values 
of the lightest heavy neutrino mass M' . Also plotted is the inverse decay rate W^P in flavour a, 
with K\ a = 3.5. The mass hierarchy has been fixed to M/M' = 10. The shaded region denotes 
the time at which the asymmetry is produced (see text). 

obtained when M/M' = 3. Increasing the hierarchy increases the scattering rate for a fixed 
value of the baryon asymmetry and thus the lower bound increases. Note that the shaded 
area here denotes the non-hierarchical region M < 3 M', where our results are not valid. 
It should be noted that this lower bound was obtained maximizing the asymmetry with 
respect to K\ and Ki a and the phases contributing to the CP asymmetry. We have chosen 
K\ a = 3.5, the smallest value in the strong washout regime 4(|. Moreover, we have fixed 
Kip^a = 11, so that K\ = 25.5, which maximizes the asymmetry. 1 

In our model leptogenesis can take place more than one order of magnitude below the 
usual lower bound, given by 3 x 10 9 GeV at the onset of the strong washout. Unfortunately, 
it turns out that the coefficient of the d = 6 operator leading to non-unitarity in Eq. (jSJ) 
is expected to be very suppressed, of order 77^ ~ 10~ 16 . The reason is mainly the high 
scale of M', and thus M, needed for successful leptogenesis. Moreover, the case M/M' = 3 



1 If is too close to |Yi a |, the washout in all flavours are similar, and the asymmetry is suppressed since 
E\ ~ 0. However, if \Yip\ is much bigger than |Yi Q |, then the CP asymmetry is suppressed [c.f. Eq. ([27)) ]. 
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FIG. 3: Lower bound on the lightest heavy neutrino mass, M^ in , for successful leptogenesis, vs. 
the hierarchy parameter M/M' . The washout parameter is fixed to K\ a = 3.5, and the asymmetry 
is mainly produced in flavour a. 

provides the largest non-unitarity possible, since the latter decreases when the scale of M 
increases. For instance, with M/M' = 10, the Yukawa coupling needed for leptogenesis 
is y ~ 0.07, in which case we find the lower bound M^ in ~ 10 9 GeV and the expected 
deviation from unitarity is rj a p ~ 10~ 18 . We conclude that, in the considered scenario, 
successful leptogenesis is incompatible with observable non-unitarity signals. 

VI. SUMMARY AND DISCUSSION 

We have investigated the relation between non-unitarity of the leptonic mixing matrix 
and baryogenesis via thermal leptogenesis. We have first studied how all parameters of the 
canonical type-I seesaw mechanism can, in principle, be reconstructed from the neutrino 
mass matrix and a measurement of the deviation of the effective low-energy leptonic mixing 
matrix from unitary. In the effective low-energy theory, neutrino masses and non-unitarity 
are encoded in the lepton number violating d — 5 (Weinberg) operator and in the d = 6 
operator contributing to the neutrino kinetic terms after electroweak symmetry breaking, 
respectively. 
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For the case that the mass M' of the lightest right-handed neutrino is lighter than the 
masses of the others, we show that its decay asymmetries for flavour-dependent leptogen- 
esis can be expressed in terms of two contributions, one depending on the unique d = 5 
operator generating neutrino masses and one depending on the d — 6 operator associated 
with non-unitarity. We have argued that in low-energy seesaw scenarios, where small lepton 
number violation explains the smallness of neutrino masses, the lepton number conserving 
d — 6 operator contribution, linked to non-unitarity, generically dominates over the d — 5 
operator contribution which results in a strong enhancement of the flavour-dependent decay 
asymmetries without any resonance effects. We have referred to this case as non-unitarity 
driven leptogenesis. 

To calculate the produced final baryon asymmetry, however, we found that lepton flavour 
equilibrating effects directly related to non-unitarity play a crucial role and their effects have 
to be included. In the simple realization of non-unitarity driven leptogenesis considered here, 
they turn out to forbid lowering the leptogenesis scale down to the TeV scale. Nevertheless, 
lowering of the leptogenesis scale, i.e., the mass of the lightest right-handed neutrino M', 
to about 10 8 GeV is possible, which is more than one order of magnitude below the scale of 
standard thermal leptogenesis. The reduced leptogenesis scale in non-unitarity driven lepto- 
genesis can improve consistency between leptogenesis and gravitino (or similar) constraints 
in supergravity theories. On the other hand, the deviation from unitarity for the case of 
M' > 10 8 GeV is far below the experimentally accessible region. 
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